Sparse Matrices and Vectors in C-XSC
نویسندگان
چکیده
C-XSC is a C++ library for reliable scientific computing, which provides data types for dense vectors and matrices with real, complex, real interval and complex interval entries. These data types are easy to use and provide many helpful functionalities such as the ability to work with submatrices and subvectors. However, when dealing with sparse vectors, and especially with sparse matrices, these data types are inefficient. CXSC version 2.4.0 added special types for sparse vectors and matrices that take advantage of the sparsity, both for performance and for memory consumption. This paper explains the data structures for and the implementation of these new types. Many examples and some performance tests with sparse matrices from real world applications are included.
منابع مشابه
An Accurate an Efficient Selfverifying Solver for Systems with Banded Coefficient Matrix
In this paper we discuss a selfverifying solver for systems of linear equations Ax = b with banded matrices A and the future adaptation of the algorithms to cluster computers. We present an implementation of an algorithm to compute efficiently componentwise good enclosures for the solution of a sparse linear system on typical cluster computers. Our implementation works with point as well as int...
متن کاملSome new restart vectors for explicitly restarted Arnoldi method
The explicitly restarted Arnoldi method (ERAM) can be used to find some eigenvalues of large and sparse matrices. However, it has been shown that even this method may fail to converge. In this paper, we present two new methods to accelerate the convergence of ERAM algorithm. In these methods, we apply two strategies for the updated initial vector in each restart cycles. The implementation of th...
متن کاملImproving the Performance of a Verified Linear System Solver Using Optimized Libraries and Parallel Computation
A parallel version of the self-verified method for solving linear systems was presented in [16, 15]. In this research we propose improvements aiming at a better performance. The idea is to implement an algorithm that uses technologies as MPI communication primitives associated to libraries as LAPACK, BLAS and C-XSC, aiming to provide both self-verification and speed-up at the same time. The alg...
متن کاملNon-Asymptotic Theory of Random Matrices Lecture 16: Invertibility of Gaussian Matrices and Compressible/Incompressible Vectors
We begin this lecture by asking why should an arbitrary n × n Gaussian matrix A be invertible? That is, does there exist a lower bound on the smallest singular value s n (A) = inf x∈S n−1 Ax 2 ≥ c √ n where c > 0 is an absolute constant. There are two reasons (or cases) which we will pursue in this lecture. 1. In Lecture 15 we saw that the invertibility of rectangular (i.e., non-square) Gaussia...
متن کاملSOLUTION-SET INVARIANT MATRICES AND VECTORS IN FUZZY RELATION INEQUALITIES BASED ON MAX-AGGREGATION FUNCTION COMPOSITION
Fuzzy relation inequalities based on max-F composition are discussed, where F is a binary aggregation on [0,1]. For a fixed fuzzy relation inequalities system $ A circ^{F}textbf{x}leqtextbf{b}$, we characterize all matrices $ A^{'} $ For which the solution set of the system $ A^{' } circ^{F}textbf{x}leqtextbf{b}$ is the same as the original solution set. Similarly, for a fixed matrix $ A $, the...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- Reliable Computing
دوره 14 شماره
صفحات -
تاریخ انتشار 2010